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An alternative solution to the 7-ray Gradient problem 
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The Fermi-LAT collaboration recently confirmed EGRET finding of a discrepancy between the observed longi- 
tudinal profile of 7-ray diffuse emission from the Galaxy and that computed with GALPROP assuming that cosmic 
rays are produced by Galactic supernova remnants. The accurate Fermi-LAT measurements make this anomaly 
hardly explainable in terms of conventional diffusion schemes. Here we use DRAGON numerical diffusion code to 
implement a physically motivated scenario in which the diffusion coefficient is spatially correlated to the source 
density. We show that under those conditions we are able to reproduce the observed fiat emissivity profile in 
the outer Galaxy with no need to change the source term, the diffusion halo height, or the CO-H2 conversion 
factor (Xco) with respect to their preferred values/distributions. We also show that our models are compatible 
with gamma-ray longitudinal profiles measured by Fermi-LAT, and still provide a satisfactory fit of all observed 
secondary-to-primary ratios, such as B/C and antiprotons/protons. 



1. Introduction 

It has been known since the EGRET era that, if 
one computes the cosmic ray (CR) Galactocentric ra- 
dial distribution adopting a source function deduced 
from pulsar or supernova remnant (SNR) catalogues, 
the result appears much steeper than the profile in- 
ferred from the 7-ray diffuse emission along the Galac- 
tic plane: the latter appears flatter, with a high con- 
tribution from large Galactic radii. This discrepancy 
is known as the j-ray gradient problem. A sharp rise 
of the conversion factor between CO emissivity and 
H2 density (the so called Xco) with the Galactocen- 
tric radius was invo ked at the time to fix the problem 
[strong et al. I l2004j : a larger gas density at large radii 
compensates for the decreasing CR population and is 
able to explain the 7-ray fiux detected at high Galactic 
longitudes. 

Fe rmi-LAT confi rmed the existence of such a prob- 
lem [Ackermann et al. 2011J. Moreover, the high spa- 
tial resolution of the LAT permitted to disentan- 
gle the emission coming from the interaction of CRs 
with the molecular gas (whose modelling is strongly 
affected by the uncertainty on the Xco) from the 
emission originated by the interaction of the Galac- 
tic CRs with the atomic gas (whose density is bet- 
ter known from its 21 cm radio emission). An anal- 
ysis based on 7-ray maps of the third Galactic quad- 
rant [Ackermann et all 12011 pointed out that the 7- 



ray emissivity from neutral gas (tracing the actual 
CR density) is indeed flatter than the predicted one 
confirming the gradient problem in dependently of the 
Xco • This result led the authors of [ Ackermann et ahl 



I2OII I to look for alternative explanations of the prob- 
lem, e.g. invoking a thick CR diffusion halo or a source 
term that becomes flatter at large radii. Both solu- 
tions, however, do not appear completely satisfactory: 
a thick halo is disfavoured both from ^"Be/^Be and 
synchrotron data; a smooth source distribution is in 
contrast with SNR catalogues. 

Here we consider a different interpretation based on 
relaxing the approximation of isotropic and spatially 
uniform diffusion. 



2. Inhomogeneous and anisotropic 
diffusion 

Nearly all CR diffusion models presented in the lit- 
erature adopt an isotropic and spatially uniform dif- 
fusion coefficient throughout all the Galaxy. This is 
the case, for example, of GALPROP numeric al pack- 
age o n which the predictions of ^Ackermann et al.l 
I2OII are based. It is reasonable, however, to ex- 
pect that CR diffusion is not isotropic in the Galaxy. 
This c ould be the consequence ei ther of Galactic 
winds [Gebauer and de Boeil l2009l | or just of the 
anisotropy of the regular component of the Galac- 
tic magnetic field which is oriented almost azimu- 
tally along the Galactic plane. The former possi- 
bility was suggested as a possible solution of the 7- 
ray gra dient problem originated by EGRET obser- 
vations [Breitschwerdt et al. I l2002{ . Here we consider 
the latter option and extend also to the recent Fermi- 
LAT data the arguments we developed in [Evoli et 
al. 2008] to interpret earlier EGRET measurements. 



eConf C110509 



2 



2011 Fermi Symposium, Roma., May. 9-12 



Our approach is based on the consideration that, 
for geometrical reasons, CRs should escape from the 
Galaxy almost perpendicularly to the Galactic plane: 
their density, therefore, should be determined by the 
perpendicular component of the diffusion coefficient 
D^. We know - both from quasi linear diffusion 
theory and frorn mor e realistic numerical simulations 
iDeMarco et al. I l2007f - that D± should increase with 
increasing strength of the Galactic magnetic field tur- 
bulent component; from a physical point of view, such 
behaviour can be understood in terms of magnetic 
field line random walk becoming stronger when the 
turbulence strength increases. As a consequence, the 
regions where CR injection is more intense should also 
be those characterized by a stronger MHD turbulence 
and hence a faster CR escape along the z axis: this 
should smooth the CR gradient, and hence the 7-ray 
profile, in a rather natural way. 

In the next section we will show this effect by means 
of dedicated numerical computations. 




3. Our method and results 




Figure 2: Two gamma-ray longitudinal profiles along the 
Galactic plane computed with DRAGON and GammaSky and 
compared to preliminary Fermi- L AT data extracted from the 
talk by A.W. Strong at the Workshop on Indirect Dark Matter 
Searches, DESY, Hamburg, June 2011 

(tttp : //www .mpe .mpg. de/~aws/talks/ ). Data are integrated 
over the latitude interval —5° < b < +5° and in energy 
between IIO4 and 144^ MeV. Red line: IC. Green line: 
Bremsstrahlung. Blue line: tt" decay. Purple line: 
contribution from unresolved sources. Grey line: + IC -f- 
Bremsstrahlung. Black line: total. Panel a) The profile is 
computed with no radial dependence of diffusion coefficient. 
Panel b) Here the diffusion coefficient follows the source 
term: D oc Q'^ , with t = 0.8. 




Figure 1: Two different CR proton distribution maps in 
arbitrary units computed with DRAGON at 10 CeV are shown as 
functions of the Galactic cylindrical coordinates R and z. 
Panel a) The proton distribution is computed with no radial 
dependence of diffusion coefficient. Panel b) Here the 
diffusion coefficient is correlated to the source term: D oc Q'^ , 
with T = 0.8. The model shows a significant flattening in the 
CR profile along R. The normalization is fixed at -Rsun = 8.5 
kpc in both cases; notice how the maximum proton density is 
reduced by a factor ~ 2 in the second panel. 



In this section we use DRAGON numerical diffusion 
package^ to solve the diffusion equation in the pres- 
ence of a diffusion coefficient spatially correlated to 
the CR source term Q{r, z). We perform our analysis 
using a Plain Diffusion (PD) setup with no convec- 
tion and no reacceleration in order to better highlight 
the effects of inhomogeneous diffusion; similar results 
may be obtained with different choices of the diffu- 
sion parameters, and a more detailed study on the 
effects of another setup will be performed in a forth- 
coming paper. The CR propagation model adopted 
he re is basically the same a s the PD model described 



Di Bernardo et aT]|201l| ; the astrophysical param- 



eters (in particular the source term, gas distribution 
and Xco) are also the same used in that analysis. 
Only the normalization of the proton injection spec- 



^The DRAGON code for cosmic-ray transport 
and diffuse emission production is available online at 
|http : //www . desy ■ de/-macclone/DRAGOH/ | 
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Emissivity Gradient 




Figure 3: Here the effect of the parameter r defined by 
eauation\l\is explored. Dotted line: no radial dependence of 
diffusion coefficient (t = 0). Dot-dashed line: t = 0.2. 
Dashed line: r = 0.5. Solid lines: r = 0.7 - 0.8 - 1.0. The 
values corresponding to the solid lines within the grey band 
match the observed gradient. 



trum is slightly tuned to match the recently released 
proton spectrum measured by the PAMELA collabo- 
ration 



Adriani et al.|[20Tl |. The model is also com- 
patible with most other CR data sets. Only a little 
excess in the antiproton spectrum must be pointed out 
which, however, is still compatible with data if astro- 
physical and particle physics uncertainties are taken 
into account. As we mentioned, the CR distribution is 
computed with DRAGON: this numerical package is suit- 
able for our purpose since, differently from GALPROP, 
it implements the possibility to vary the diffusion co- 
efficient through the Galaxy. The CR distributions 
is then used as an input to compute the 7-ray longi- 
tude profile along the Galactic plane; the 7-ray map 
is evaluated with a separate package called GammaSky. 

The result of a combined DRAGON and GammaSky 
computation in the case of a uniform diffusion co- 
efficient and a PD setup is shown in Fig. [5] (panel 
a). It is clear from that plot that the predicted lon- 
gitude profile is too steep compared to the observa- 
tions: in the Galactic center region the model predic- 
tion overshoots the data and in the anti-center region 
the model is lower than the observations by several a. 
Tuning the Xco(^) could help in principle; assuming 
a lower value of this parameter in the bulge and a high 
value at large R could smooth the 7-r ay profile (as 
done in several previous works such as [Strong et al.l 
|2004 |). Unfortunately, as pointed out in the introduc- 
tion, the gradient problem is present especially in the 
emissivity profile, and this quantity is independent of 
the molecular gas: it only traces the actual CR distri- 
bution^ . 

So we apply our previous considerations and adopt 



a diffusion coefficient correlated to the radial depen- 
dence of the source term Q{R) by the following ex- 
pression: 

D{R) oc QiRY (1) 

This is the p arametrization we already used in 
[Evoli et all [2C)0&] to interpret EGRET data. The 
parameter r is tuned against data. In Fig. [3] we 
show the emissivity profile for different values of r in 
the range [0^1]. It is evident from that figure that 
an increasing value of r yields a much smoother be- 
haviour of the emissivity as function of R. Values in 
the ra nge [0.7 1] allow a good match of Fermi-L AT 
data f [Ackermann et al.ll201l| . jAbdo et al.ll2010l |). 

With this result at hand, we considered a modified 
version of the Plain Diffusion CR propagation setup 
with D{R) = Q'^ and t = 0.8. The smoothing in the 
CR distribution corresponding to such a value of t 
is shown in Fig. [T] As shown in Fig. [5] (panel b), 
the 7 ray longitude profile along the Galactic plane is 
nicely reproduced with no tuning at all of the Xco- 
It is remarkable that a simple CR propagation setup, 
with only the addition of the radial dependence of D 
and no ad hoc tuning, permits to reproduce the 7- 
ray profile with such accuracy. Noticeably, the mod- 
ified model is still compatible with most relevant CR 
data set, most importantly the B/C. Furthermore, we 
checked that the 7-ray spectrum measured by Fermi- 
LAT along the Galactic plane is also correctly repro- 
duced under those conditions (see Fig. [1]). 




^The emissivity is the number of 7 photons emitted by each 
gas atom per unit time and unit energy 



Figure 4: The gamma-ray spectrum corresponding to the 
plain diffusion model with varying diffusion coefficient 
described in the text (r = 0.8j. The spectrum was computed 
with DRAGON and GaimnaSky. The data points measured by 
Fermi-LAT are taken from the same reference as Fig. [2| 



4. Conclusions 

In this paper we presented an alternative solution to 
the well known 7-ray gradient problem. Our approach 
is based on the physically motivated hypothesis that 
the CR diffusion coefficient is spatially correlated 
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to the source density: regions in which star, hence 
SNR, formation is stronger are expected to present a 
stronger turbulence level and therefore a larger value 
of the perpendicular diffusion coefficient. This effect 
favours CR escape from most active regions helping to 
smooth their density through the Galaxy hence also 
the 7-ray gradient. We used DRAGON package to im- 
plement this scenario and to check that CR data are 
still reproduced under those conditions. In spite of be- 
ing purely phenomenological (as a self-consistent the- 
ory/computation of non-linear CR - MHD turbulence 
interaction in the Galaxy is far from being developed) 
our approach provides a remarkably good description 
of the spectrum and longitude distribution of the dif- 
fuse 7-ray emission measured by the Fermi-LAT col- 
laboration. 
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